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Abstract 


Recent  theoretical  and  observational  studies  have  indicated 
that  neutral  wind  motions  in  the  earth's  ionospheric  E-region 
play  an  important  role  in  the  formation  and  vertical  motion  of 
thin,  concentrated  layers  of  metallic  ions.  These  metallic  ions 
are  commonly  called  sporadic-E  or  Es.  The  formation  of  sporadic-E 
layers  can  be  explained  by  considering  a steady-state  equilibrium 
between  neutral  wind  motions  and  the  resultant  collisional  and 
geomagnetic  forces  experienced  by  the  metallic  ions  in  the  neutral 
gas. 

The  modal  oscillations  of  atmospheric  tidal  winds  generate 
distinct  patterns  of  ion  convergence  which  control  the  formation 
and  movement  of  Es  layers.  The  interaction  of  thermally-excited 
tidal  winds  and  metallic  ions  in  the  presence  of  the  earth's  mag- 
netic field  has  been  numerically  simulated.  A comparison  of  the 
simulation  predictions  with  observational  evidence  of  sporadic-E 
occurrence^  from  two  widely-  separated  iononsonde  stations  indicates 
that  tidal  winds  play  a distinct  and  predominant  role  in  governing 
the  behavior  of  mid-latitude  Es.  The  initial  results  suggest  that 
it  may  be  possible  to  predict  the  gross  synoptic  features  of  spora- 
dic-E on  a global  scale. 
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I.  INTRODUCTION 

Sporadic-E,  often  designated  E$  , is  a term  used  to  describe 
thin,  highly-concentrated,  horizontally-stratified  layers  of  ioniza- 
tion which  frequently  form  In  the  earth's  ionospheric  E-reglon  (100- 
150  km).  The  adjective  "sporadic"  was  originated  by  early  Ionospheric 
researchers  who  recognized  that  E$  behavior  is  not  controlled  by 
classical  Chapman  (1931a,  1931b)  photochemical  mechanisms. 

In  general,  sporadic-E  layers  have  a common  set  of  character- 
istics. The  vertical  thickness  of  most  E$  layers  typically  ranges 

between  100  meters  and  5 kilometers.  Peak  ionization  concentrations 

4 6 -3 

are  usually  between  10  and  10  electrons  cm  . Therefore,  the 

vertically-integrated  total  electron  content  for  sporadic-E  layers 

8 10  - 2 

spans  the  range  from  10  to  10  electrons  cm  (MacLeod,  1975).  Since 

the  normal  photochemically-produced  ionization  in  the  E-region  is  only 
3 5 -3 

on  the  order  of  10  -10  electrons  cm  , sporadic-E  layers  appear  as 
distinct  discontinuities  in  ambient  ionization  concentration  profiles. 

Several  typical  sporadic-E  electron  concentration  profiles 
are  depicted  in  Figure  1.  The  series  of  profiles  shown  in  the  figure 
illustrate  several  important  features  which  are  common  to  many 
sporadic-E  layers.  First,  sporadic-E  layers  tend  to  form  in  the  upper 
reaches  of  the  ionospheric  E-region,  often  between  150  and  175  kilo- 
meters. Subsequently,  they  descend  to  approximately  110  kilometers 
over  a period  ranging  from  6 to  12  hours.  The  descent  is  often 
accompanied  by  an  enhancement  of  the  ionization  concentration  of  the 
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layer.  In  most  cases,  the  Increase  in  the  ionization  concentration 


appears  to  be  caused  by  a vertical  compression  of  the  ionized  material 
contained  within  the  layer  rather  than  by  the  addition  of  new  ioniza- 
tion. Frequently,  more  than  one  sporadic-E  layer  is  present  in  the 
E-region  at  the  same  time  (see  Figure  1,  profile  1).  More  often  than 
not,  the  lower  E$  layer  is  a remanent  of  an  older  layer  which 
descended  from  the  upper  E-region  several  hours  earlier.  The  tendency 
for  the  pattern  of  descending  E$  layers  to  repeat  itself  in  a 
cyclical  fashion  has  earned  this  behavior  the  name  "sequential 
sporadic-E".  To  date,  the  reverse  process,  i.e.,  ascending  E$  layers, 
has  never  been  reported.  A detailed  discussion  of  sequential  sporadic- 
E,  as  well  as  the  role  played  by  atmospheric  tidal  winds  in  its  forma- 
tion and  propagation,  will  be  presented  in  the  subsequent  chapters  of 
this  report. 


Figure  1.  Night-time  sequence  of  electron  density  profiles  as 
observed  at  Wallops  Island,  Virginia,  by  Smith  (1970) 
(taken  from  Smith,  1970). 
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Systematic  observations  of  sporadic-E  Indicate  that  the 


layers  often  cover  large  geographical  areas.  E$  "clouds"  spanning 
horizontal  distances  of  more  than  500  kilometers  have  been  measured 
by  sporadlc-E  lonosonde  networks  In  Europe  (Bossy,  1972).  Other  less 
reliable  measurements  (Gersen,  1955)  suggest  that  the  horizontal  ex- 
tent of  the  Ionized  layers  may  exceed  1000  kilometers.  At  the  same 
time,  most  rocket  and  lonosonde  observations  Indicate  that  the 
ionized  material  contained  within  large  sporadic-E  layers  Is  seldom 
uniformly  distributed.  In  particular,  order  of  magnitude  variations 
in  Ionization  concentrations  of  the  layers  over  horizontal  and  verti- 
cal distances  of  less  than  1 kilometer  appear  to  be  the  rule  rather 
than  the  exception  (Bossy,  1972).  Clearly,  the  extreme  Inhomogenetles 
in  the  distribution  of  ionization  within  sporadic-E  clouds  suggest 
that  complex  microscale  mechanisms  are  operating  within  E$  layers. 
However,  it  is  not  the  intent  of  this  paper  to  examine  the  small- 
scale  ((ynamics  of  sporadic-E  layers.  Only  large,  global-scale  features 
will  be  considered.  None-the-less,  the  complex  internal  structure  of 
the  layers  will  necessarily  complicate  comparisons  and  limit  the  de- 
gree to  which  large-scale  theoretical  model  predictions  can  be  verified 
using  conventional  lonosonde  data. 

Rocket  and  ground-based  measurements  of  the  chemical  composi- 
tion of  sporadlc-E  layers  (Istomin,  1963;  Young  et.  al . , 1967); 

Goldberg  et.  al.,  1973;  Philbrlck  et.  al . , 1974;  Behnke  et.  al.,  1975) 
all  clearly  indicate  that  the  layers  are  largely  composed  of  metallic 
ions.  Although  Fe^g  , SI , and  Mg^  are  usually  the 
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predominant  ions,  other  elemental  ions,  notably  Nigg  , Cr^  » Ca^g  » 
and  Na^  » have  been  detected  with  regularity  and  considerable  cer- 
tainty. The  presence  of  metallic  elements  in  sporadic-E  layers  in 
roughly  the  same  proportions  as  the  elemental  composition  of  meteors 
suggests  that  E$  layers  are  composed  of  either  finely-divided 
meteoric  dust  or  debris  ablated  from  larger  meteors.  The  outstanding 
correlation  between  meteor  showers  and  the  intensity  and  frequency  of 
occurrence  of  sporadic-E  (Appleton  et.  a!.,  1947)  lends  strong  addi- 
tional support  to  the  importance  of  the  role  played  by  meteoric 
material  in  the  production  of  sporadic-E. 

The  presence  of  metallic  elements  and  compounds  in  sporadic-E 
accounts  for  several  unique  features  of  E$  behavior.  First, 
metallic  elements  and  metallic  oxide  compounds  have  very  low  ioniza- 
tion and  dissociation  potentials  (Swider,  1969;  Opik,  1958).  There- 
fore, meteoric  material  is  readily  ionized  through  charge  exchange, 
photoionization  or  aerodynamic  heating  processes.  Second,  once 
metallic  ions  are  formed,  they  are  lost  very  slowly.  The  primary 
kinetic  and  chemical  loss  mechanisms  include  three  body  collisions 
and  radiative  recombination.  Recombination  coefficients  for  either 
process  are  typically  of  the  order  10"^®  cm"*  sec"^  . Therefore, 
lifetimes  for  the  more  abundant  metallic  Ions  are  approximately  one 
day  in  the  lower  E-region  (Whitehead,  1970). 

The  unique  properties  of  the  ionized  metallic  elements  and 
compounds  found  in  sporadic-E  layers  create  a regime  in  which  Ions 
are  rapidly  produced  and  slowly  lost.  Consequently,  the  behavior 


4 


of  sporadic-E  layers  in  general  and  metallic  ions  in  particular  is 
largely  dominated  by  wind-driven  and  electrostatic  transport  mechanisms 
while  photochemical  processes  play  a secondary  role.  Accordingly, 
many  ionospheric  researchers  have  attempted  to  identify  the  dynamic 
atmospheric  transport  mechamisms  which  can  account  for  the  observed 
formation  and  vertical  propagation  of  sporadic-E  layers. 

Presently,  it  is  widely  recognized  that  vertical  shears  in 
horizontal  ionospheric  winds  play  an  important,  if  not  predominant,  rol 
in  the  production  of  mid-latitude  sporadic-E.  A large  theoretical 
body  of  knowledge  appropriately  termed  the  "sporadic-E  wind-shear 
theory"  has  been  developed  and  expanded  upon  by  numerous  ionospheric 
physicists,  notably,  Dungey  (1956,  1959),  Whitehead  (1961),  Smith  and 
Matsushita  (1962),  MacLeod  (1964,  1966,  1968,  1969),  MacLeod  et.  al . , 
(1973),  1975),  Keneshea  et.  al . , (1970).  Furthermore,  the  validity 
and  importance  of  wind-driven  ion  convergence  mechanisms  have  been 
repeatedly  confirmed  by  comprehensive  theoretical  and  observational 
studies  (MacLeod  1966,  1975).  The  success  of  wind-driven  sporadic-E 
theories  in  explaining  the  formation  of  many  E$  layers  should  not 
rule  out  the  possibility  that  other,  more  complex,  sporadic-E 
mechanisms  are  not  concurrently  functioning.  Indeed,  the  inability  of 
the  wind-shear  theory  to  explain  all  the  occurrences  and  features  of 
mid-latitude  sporadic-E  suggests  that  all  of  the  important  formation 
mechanisms  have  not  yet  been  discovered.  However,  speculation  regard- 
ing the  importance  of  sporadic-E  mechanisms  which  are  driven  by 
forces  other  than  neutral  wind  motions  and  their  associated  dynamo 


electric  fields  is  beyond  the  intended  scope  of  this  paper.  Further- 
more, in  order  to  eliminate  the  additional  E-region  complications 
brought  about  by  energetic  particle  precipitation  at  high  latitudes 
and  the  equatorial  electrojet  at  low  latitudes,  the  subject  matter  in 
this  thesis  will  be  further  confined  to  the  behavior  of  mid-latitude 
sporadic-E  only. 

The  specific  purpose  of  this  thesis  is  to  examine  in  detail  the 
role  played  by  atmospheric  tidal  winds  and  electrostatic  fields  in 
the  formation  and  vertical  propagation  of  ionospheric  sporadic-E  layers. 
Although  numerically-oriented  sporadic-E  models  have  been  developed 
by  other  researchers  (MacLeod,  1975)  in  order  to  compare  rocket-derived 
wind  measurements  with  observed  occurrences  of  sporadic-E,  there  are 
apparently  no  published  accounts  of  attempts  to  numerically  model 
and  interpret  the  behavior  of  sporadic-E  layers  produced  by  atmospheric 
tidal  winds.  Therefore,  it  may  be  concluded  that  some  of  the  results 
presented  in  this  paper  are  somewhat  new  and  potentially  useful  in 
predicting  periodic  features  of  mid- latitude  sporadic-E  layers. 

For  the  reader's  benefit,  this  thesis  has  been  divided  into 
several  functionally- independent  chapters.  The  intended  objectives 
of  each  chapter  can  be  briefly  summarized  as  follows: 

Chapter  I:  Introduction. 

a)  Define  sporadic-E  and  present  a brief  outline 
of  its  interesting  features. 

b)  Highlight  the  importance  of  metallic  ions  and 
ionospheric  neutral  winds  in  the  formation  of 
Eg  layers. 
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Chapter  II:  Analysis  of  wind-driven  ion  motion  and 
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a) 


b) 

Chapter  III: 


a) 


b) 


Chapter  IV. 


convergence. 

Develop  a mathematical  analysis  of  the  impor- 
tant wind-driven  and  electrostatic  forcing 
mechanisms  which  are  capable  of  producing 
sporadic-E  layers. 

Offer  a detailed  interpretation  of  the  physical 
implications  of  the  analysis. 

The  theoretical  role  of  tidal  winds  in  the  pro- 
duction and  vertical  propagation  of  mid-latitude 
sporadic-E. 

Discuss  how  the  mathematical  analysis  presented 
in  Chapter  II  has  been  combined  with  numeri- 
cal models  of  the  earth's  magnetic  field, 
upper  air  density  and  thermally-excited 
atmospheric  tides  for  the  purpose  of  analyzing 
the  influences  of  both  tidal  winds  and  electro- 
static fields  on  the  formation  and  behavior 
of  sporadic-E. 

Illustrate  and  discuss  the  theoretical  spatial 
and  temporal  sporadic-E  patterns  which  should 
be  produced  by  the  dominant  (equinox) 
atmospheric  tidal  wind  modes. 

Comparison  of  theoretical  results  with  sporadic-E 
data. 
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a)  Compare  the  theoretically-derived,  ttdally- 
driven  sporadic-E  patterns  from  Chapter  III 
with  high-resolution  sporadic-E  observations 
from  MacDougall  (1974a, b)  and  others. 

No  detailed  review  of  sporadic-E  research  is  attempted  in 
this  paper.  However,  the  reader  may  find  the  following  references 
especially  useful  in  this  regard: 

1)  Ionospheric  Sporadic-E,  Smith  and  Matsushita, 
MacMillan  Co.,  N.Y.,  1962. 

2)  Radio  Science,  v.  T_t  n.  3,  March  1972,  Special  Issue 
on  Sporadic-E. 

3)  Radio  Science,  v.  1£,  n.3,  March,  1975,  Special  Issue 
on  Sporadic-E. 

4)  Whitehead,  J.D.  Production  and  prediction  of 
sporadic-E,  Reviews  of  Geophysics  and  Space  Physics, 
8:65,  February  1970. 
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II.  ANALYSIS  OF  WIND-DRIVEN  ION  MOTION  AND  CONVERGENCE 


As  discussed  in  the  Introduction  to  this  paper,  upper  atmos- 
pheric neutral  winds  and  long-lived  metallic  ions  have  been  observed 
to  play  an  important  part  in  the  formation  of  ionospheric  sporadic-E 
layers.  The  purpose  of  this  chapter  is  to  mathematically  analyze  the 
manner  in  which  atmospheric  winds,  working  in  conjunction  with  magnetic 
and  collisional  forces,  can  compress  metallic  ions  into  thin, 
horizontally-stratified,  sporadic-E  layers. 

In  general,  the  earth's  ionospheric  E-region  can  be  considered 
a weakly-ionized  and  magnetized,  collisional  plasma.  In  the  central 
portion  of  the  E- region  (=  125  km)  , the  ion  gyrofrequency,  tUj  , 
for  moderately  heavy  ions  (=  35  a.m.u.)  can  be  shown  to  be  approxi- 
mately equal  to  the  ion-neutral  collision  frequency,  v-n  . (See 
Figure  2,  profile  1.)  In  the  lower  reaches  of  the  E-region  (*  100  km), 
vin  a factor  of  10  larger  than  caj  • In  the  upper  ex- 

tents of  the  ionospheric  E-region  (*  150  km)  just  the  opposite  is  true: 
v1>n  is  a factor  of  10  smaller  than  . Therefore,  for  the  most 
part,  the  ions  can  be  considered  to  be  collisionally  coupled  to  neutral 
(uncharged)  atmospheric  constituents  of  the  E-region.  On  the  other 
hand,  since  the  electron  gyrofrequency  Is  nearly  1800  times  larger 
than  the  ion  gyrofrequency , the  free  electrons  in  the  ionospheric  E- 
region  are  largely  decoupled  from  the  neutrals.  Therefore,  forced, 
bulk  motions  of  the  neutral  E-region  constituents  are  quickly 


I 
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transmitted  to  the  ions,  and  not  the  electrons,  through  collisions. 


Accordingly,  the  wind  driven  motion  of  electron-ion  pairs  in  the 
E-region  is  largely  governed  by  the  vector  equation  of  motion  for  the 
ions: 


min1  ar  = m1ni  v1n  ^ -V  + Vi  + + mini  a 


where  m^ : 

V 

vin 

U: 


V 

1- 


ion  mass 

ion  concentration 
ion-neutral  collision  frequency 
velocity  of  the  neutrals 
velocity  of  the  ions 
ion  charge 
gravity  vector 

ion  partial  pressure  gradient  (V^KI^ )) 


The  objective  at  this  point  is  to  introduce  reasonable  assump- 
tions, considerations  and  observations  to  simplify  the  process  of 
solving  the  equation  of  motion  for  the  vector  ion  velocity,  V,.  . 
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First,  it  appears  reasonable  to  assume  that  ion  partial  pres- 
sure gradients  In  the  Ionosphere  are  trivially  small  In  comparison  to 
the  other  forces  In  the  equation  of  motion.  This  Is  undoubtably  the 


) 


case  In  classical,  non-metalllc,  Chapman  Ionospheric  layers  where 
photochemical  processes  clearly  predominate  over  transport  mechanisms. 
The  behavior  of  metallic  Ion  layers  is  considerably  different.  The 
highly-stable  chemical  characteristics  of  metallic  ions  permit  wind- 
driven  transport  to  predominate  over  the  subsidiary  role  of  photo- 
chemistry. Therefore,  It  is  possible  that  Ion  partial  pressure 
gradients  grow  to  significant  proportions  in  extremely  thin 
(<  100  meters)  sporadic-E  layers  where  strong  vertical  congressional 
forcing  is  evident. 

Unfortunately,  the  inclusion  of  the  Ion  partial  pressure 
gradient  term  in  the  equation  of  motion  necessitates  the  concurrent 
consideration  of  the  ion  continuity  equation  since: 

3n. 

yjr-  = Q.j  - L.j  + V*(DVn.j  - n^)  (2) 


where : 


0: 


ni 


ion  source  function 

ion  loss  function 

ambi polar  diffusion  coefficient 

as  previously  defined 
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It  is  essential  that  the  ion  partial  pressure  gradient  term 
be  included  in  any  numerical  model  written  for  the  expressed  purpose 
of  accurately  determining  absolute  metallic  ion  concentrations  in 
sporadic-E  layers.  However,  since  the  objective  of  this  thesis  is  to 
study  only  transport  mechanisms  and  not  to  determine  absolute  metallic 
ion  concentrations,  the  ion  partial  pressure  gradient  term  will  be 
omitted  without  further  consideration.  This  assumption  effectively 
decouples  the  equation  of  motion  from  the  ion  continuity  equation, 
eliminates  the  variable  n^  (ion  concentration),  and  significantly 
simplifies  the  method  of  solving  the  equation  of  motion  for  V_^  , 
the  ion  velocity. 

Next,  the  force  of  gravity  in  the  equation  of  motion,  when 
compared  to  the  collisional  and  Lorentz  forces  generated  by 
atmospheric  tidal  winds  (|U|  = 10  to  100  m/s),  can  be  considered 
insignificant.  For  example,  collisional  accelerations  (force  per 
unit  mass)  generated  by  a 10  m/s  difference  between  the  neutral  and 
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ion  velocities  (|U  - V^|)  represents  a 1000  m/s  acceleration  at 
120  kilometers  where  v1n  ~ 100  sec”^  . Therefore,  the  omission  of 
gravitational  forcing  from  the  equation  of  motion  will  result  in 
errors  on  the  order  of  only  lit.  Consequently,  the  force  of  gravity 
can  be  Ignored  in  the  context  of  the  problem  under  consideration. 


Finally,  the  ion  equation  of  motion  can  be  further  simplified 
by  considering  only  steady-state  equilibrium  solutions.  However,  in 
order  to  justify  the  use  of  steady-state  solutions,  it  is  necessary 
to  consider  the  time  required  for  mechanically-forced  ion  motions 
to  reach  a steady  state.  This  problem  has  been  examined  in  detail  by 
Piddington  (1954).  For  a mechanically-excited,  collisional,  magnetic, 
gas-ion  system,  Piddington  has  shown  hat  the  time  required  for  the 
oscillatory  ion  velocity  components  to  damp  one  e-fold  is  of  the  order 


(3) 


where  H:  magnetic  field  intensity 

-20 

: Hall  conductivity  (*  10  e.m.u.  @ 120  km) 

p:  density  of  the  neutrals  (*  5 x 10 gm/cc  @ 120  km) 

For  the  middle  of  the  ionospheric  E-region  (~  120  km),  a 
is  much  less  than  one  inverse  second.  Therefore,  the  steady-state 
ion  velocity  assumption  is  an  excellent  approximation  for  the  diurnal 
and  semidiurnal  mechanical  forcing  motions  generated  by  the  predominant 
atmospheric  tidal  wind  modes. 

By  combining  the  simplifying  approximations  discussed  in  the 
previous  paragraphs,  the  ion  equation  of  motion.  Equation  (1),  can  be 
simplified  to: 
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/ x B \ 

"1  % <Ji  - ij)  * (l  *-TT-j=  0 


This  simplified  equation  of  motion  can  be  rewritten  as; 


0i  (U  - V.)  + £ + (V.  x b)  = 0 , 


where:  0^  = 


v1n/“1 : 


w: 


ratio  of  ion-neutral  collision  frequency 

to  ion  gyrofrequency 

unit  vector  parallel  to  the  magnetic 

field 

reduced  electric  field  (dimensions  of 
velocity) 


Basically,  this  is  a statement  that  a balance  exists  between 
the  ion-neutral  colli  si onal  force  and  the  Lorentz  force. 

Through  a series  of  relatively  straight-forward  vector  manipu- 
lations, Equation  (5)  can  be  solved  for  V^-  • The  mathematical  steps 
of  the  derivation  are  presented  in  detail  in  Appendix  1.  The  essential 
results  are: 


11  •&][■** 

•bit- 


+ 0^  (U  x b)  + (U*b)  b 


I + Ixb  + [-  (1-b)  b 


Recall  that  in  Equation  (6)  is  the  steady-state  velocity 
of  an  ion  in  a mechanically-  and  electrostatically- forced,  collisional, 
magnetic,  weakly-ionized  plasma.  In  the  context  of  ionospheric 
sporadic-E,  one  is  interested  in  determining  the  circumstances  in  which 
the  ion  velocity  converges , since  ion  velocity  convergence  corresponds 
to  increasing  ionization  concentrations.  This  convergence  (negative 
divergence)  can  be  directly  computed  from  Equation  (6). 

The  complete  details  of  the  computations  of  the  ion  velocity 
divergence/convergence  from  Equation  (6)  are  presented  in  Appendices 
2 and  3.  In  light  of  the  complexity  of  the  results,  the  solution  will 
be  presented  and  discussed  in  two  separate  parts.  The  first  part  will 
cover  wind-driven  ion  velocity  convergence  (negative  divergence  of  the 
first  three  terms  on  the  right  hand  side  of  equation  (6)).  The  second 
part  will  consider  electrostatically-driven  ion  velocity  convergence 
(negative  divergence  of  the  last  three  terms  on  the  right  hand  side 
of  Equation  (6)). 


Mind-driven  ion  velocity  convergence: 

From  Appendix  2,  the  results  of  the  wind-driven  ion  velocity 
convergence  computations  are: 
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3? 

n 


1+6^ 


(V-U)  - V 


6i 


1+8? 


(U  x b) 


'*i‘ 

■ - 

b- (V  x U ) - 

1 

b-(b-v)U  - (U-b)  V 

1 

— 7 

1+6? 

.1+4 

(7) 
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Equation  (7)  implies  that  six  separate  wind-driven  mechanisms 
are  capable  of  producing  ionospheric  sporadic-E.  Note  that  the  inter- 
pretation of  this  equation  is  considerably  simplified  if  the  following 


assumptions  are  made: 

1)  Vertical  gradients  are  much  larger  than  horizontal 
gradients  (i.e.,  V = |j  ) . 

2)  Bulk  neutral  wind  motions  are  largely  horizontal 
(i.e.,  Uz  * 0). 

If  these  assumptions  are  made,  it  can  be  deduced  from  Equa- 
tion (7)  that  vertical  ion  velocity  convergence  (sporadic-E)  should 
occur  under  the  following  circumstances: 

1)  The  cross  product  of  the  neutral  wind  vector  with  the 
earth's  magnetic  field  vector  has  a vertical  component  (with  the 
appropriate  sign).  (Term  3,  Equation  (7)). 

2)  The  neutral  wind  component  perpendicular  to  the 
horizontal  projection  of  the  earth's  magnetic  field  vector  has  a 
vertical  shear  (of  the  appropriate  sign).  (Term  4,  Equation  (7)). 

3)  The  neutral  wind  vector  component  parallel  to  the 
horizontal  projection  of  the  earth's  magnetic  field  has  a vertical 
shear  (of  the  appropriate  sign).  (Term  5,  Equation  (7)). 

4)  Tne  neutral  wind  has  a vector  component  parallel  to 
the  horizontal  projection  of  the  earth's  magnetic  field.  (Term  6, 
Equation  (7)). 

Under  the  imposed  assumptions  (V  =—■  , IJ2  = 0)  , the  first 
and  second  terms  in  Equation  (6)  are  incapable  of  generating  vertical 
ion  velocity  convergence. 
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It  is  interesting  to  note  that  the  four  wind-driven  E$  pro- 
duction mechanisms  can  be  divided  into  two  classes.  The  first  class 
consists  of  sporadic-E  production  mechanisms  which  depend  on  vertical 
shears  in  the  horizontal  neutral  winds  (terms  4 and  5,  Equation  (7)). 
Henceforth,  these  terms  will  be  referred  to  as  the  "wind-shear  Eg 
mechanisms".  The  second  class  contains  the  vertical  ion  velocity 
convergence  mechanisms  which  function  because  of  vertical  gradients 
in  3^  » the  ratio  of  the  ion-neutral  collision  frequency  to  the  ion 
gyrofrequency  (terms  3 and  6,  Equation  (7)).  These  two  terms  will 
be  designated  the  "non-shear  E$  mechanisms". 

It  is  important  to  note  that  two  of  the  E$  mechanisms 
operate  only  when  the  neutral  wind  has  a component  perpendicular  to 
the  horizontal  projection  of  the  earth  s magnetic  field  (terms  3 and 
4,  Equation  (7)).  In  contrast,  the  other  two  mechanisms  (terms  5 and 
6,  Equation  (7))  function  only  when  the  neutral  wind  has  a component 
parallel  to  the  horizontal  projection  of  the  earth's  magnetic  field. 

The  geometrical  atmospheric  factors  for  the  various  terms  in 
Equation  (7)  have  been  numerically  computed  using  Jacchia's  (1971) 
neutral  density  model  and  Chapman's  (1956)  ion-neutral  collision  fre- 
quency model.  The  computations  were  based  on  a magnetic  field  Intensity 
of  .30  Gauss  and  an  ion  mass  of  24  a.m.u.  (Mg^).  The  results  of  the 
computations  are  depicted  in  Figure  2.  Profiles  3,  4,  5,  and  6 
Illustrate  how  the  geometrical  atmospheric  factors  for  terms  3,  4,  5, 
and  6 (respectively)  in  Equation  (7)  vary  with  altitude. 
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Figure  2 Geometrical  atmospheric  factors  for 
selected  ion  convergence  mechanisms. 
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A scale  analysis  of  the  important  wind  (and  electrostatic) 
ion  convergence  mechanisms  has  been  performed  using  the  coefficients 
from  Figure  2 in  conjunction  with  representative  ionospheric  (tidal) 
winds.  The  detailed  results  are  presented  in  Appendix  5.  The  results 
confirm  the  fact  that  atmospheric  geometric  effects  (P^)  play  an 
important  role  in  establishing  the  relative  importance  of  the  various 
ion  convergence  mechamisms.  In  particular,  the  following  generaliza- 
tions can  be  made: 

a)  Ion  convergence  is  dominated  by  E-W  wind  shear 
mechanisms  in  the  lower  E-region  (120-125  km)  while  the  effects  of 
N-S  wind  shears  predominate  in  the  upper  E-region  (140-150  km). 

b)  Wind  shear  ion  convergence  mechanisms  typically 
dominate  the  non-shear  mechanisms  by  factors  ranging  from  3 to  5. 

Finally,  it  should  be  recognized  that  the  vertical  Ion  velocity 
convergence  mechanisms  under  discussion  are  functions  of  the  horizontal 
component  of  the  earth's  magnetic  field.  In  particular,  the  effective- 
ness of  all  of  the  wind-driven  convergence  mechanisms  (shear  and  non- 
shear alike)  diminishes  markedly  in  the  magnetic  polar  regions  and 
increases  where  there  are  large  horizontal  magnetic  anomalies  in  the 
earth's  magnetic  field  (e.g.,  the  Far  East).  In  general,  observations 
support  these  theoretical  points  to  a high  degree  (Smith  et.  a!.,  1960). 
A more  detailed  analysis  of  the  role  played  by  magnetic  field  variations 
in  sporadic-E  formation  will  be  considered  in  the  next  chapter. 
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Electrostatically-driven  ion  velocity  convergence: 

The  ion  velocity  convergence  (negative  divergence)  produced 
by  electrostatic  fields  can  be  computed  by  taking  the  convergence  of 
the  last  three  terms  of  the  right-hand  side  of  Equation  (6).  The 
details  of  the  derivation  are  presented  In  Appendix  3.  The  essential 
results  are: 
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Equation  (8)  implies  that  six  separate  electrostatically- 
driven  mechanisms  are  capable  of  producing  ionospheric  sporadic-E. 
The  interpretation  of  Equation  (8)  is  simplified  if  the  following 
assumptions  are  made: 
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1)  Ionospheric  plasmas  are  neutral,  i.e.,  net  charge  Is 
zero  (V*£  = 0). 

2)  Time-dependent,  magnetic  field  variations  are  small. 

V X £.  a 0 . 

A 

3)  Field-aligned  electric  fields  are  small  (£/b  = 0). 

4)  Vertical  gradients  are  much  larger  than  horizontal 
gradients  (V  = |j). 

If  these  assumptions  are  made.  Equation  (8)  implies  that 
vertical  ion  velocity  convergence  (sporadic-E)  should  occur  under  the 
following  circumstances: 

1)  £ has  a vertical  component  (pointing  in  the  proper 
direction).  (Term  1,  Equation  (8)). 

2)  The  cross-product  between  £ and  the  earth's 
magnetic  field  vector  has  a vertical  component 
(pointing  in  the  proper  direction).  (Term  3, 

Equation  (8)).  This  is  the  familiar  ExB  drift. 

Under  the  imposed  assumptions,  the  remaining  terms  in  Equation 
(8)  (terms  2,  4,  5,  and  6)  are  incapable  of  generating  vertical  ion 
velocity  convergences. 

A scale  analysis  of  the  important  electrostatic  (and  wind- 
driven)  ion  convergence  mechanisms  is  presented  in  Appendix  5.  The 
results  show  that  electrostatically-driven  ion  convergence  is  typically 
an  order  of  magnitude  smaller  than  wind-driven  convergence. 


Therefore,  It  appears  that  electrostatic  fields  play  a subordinate 
role  to  the  wind-driven  convergence  mechanisms.  This  conclusion  is 
supported  by  the  work  of  MacLeod  (1964,  1975)  and  others.  Therefore, 
the  remainder  of  this  thesis  will  concentrate  on  the  major  role 
played  by  neutral  wind  motions  in  the  production  of  sporadic-E. 
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In  this  chapter,  the  velocity  of  an  ion  in  a wind-driven, 
collisional,  magnetized,  weakly- ionized  plasma  was  computed.  In  turn, 
this  result  was  used  to  compute  the  ion's  vertical  velocity  convergence. 
Subsequent  analysis  showed  that  four  wind-driven  ion  convergence 
mechanisms  are  capable  of  producing  sufficient  vertical  ion  convergence 
to  produce  sporadic-E  layers.  Additionally,  it  was  shown  that  two 
electrostatically-driven  convergence  mechanisms  can  also  contribute 
to  ion  convergence.  These  mathematical  results  will  be  used  in  the 
next  chapter  to  compute  the  sporadic-E  patterns  which  should  be 
produced  by  atmospheric  tidal  winds. 
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III.  THE  THEORETICAL  ROLE  OF  ATMOSPHERIC  TIDAL  WINDS  IN  THE 
PRODUCTION  AND  VERTICAL  PROPAGATION  OF  MID-LADITUDE 
SPORADIC-E 


The  purpose  of  the  previous  chapter  was  to  derive  the  equations 
which  describe  the  motions  and  convergence  of  metallic  Ions  In  the 
Ionospheric  E-reglon.  This  chapter  will  describe  how  these  equations 
have  been  incorporated  Into  a computer  model  of  the  earth's  Ionosphere 
for  the  purpose  of  studying  the  spatial  and  temporal  patterns  of 
sporadic-E  produced  by  various  atmospheric  tidal  wind  modes. 

Since  sporadic-E  Is  produced  largely  by  the  vertical  conver- 
gence of  metallic  ions,  the  mathematical  modelling  of  E$  layers  is 
greatly  simplified  if  the  analysis  is  restricted  to  motions  and 
convergence  along  the  vertical  coordinate  only.  Therefore,  the 
numerical  sporadic-E  model  which  will  be  presented  in  this  chapter 
is  purely  one-dimensional.  For  the  one-dimensional  case,  the  wind- 
driven,  vertical  component  of  the  ion  velocity  equation  can  be  written 
as: 
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Next,  the  vertical  ion  velocity  convergence  can  be  readily 
computed  from  Equation  (9): 
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In  order  to  numerically  simulate  the  behavior  of  sporadic-E, 
mathematical  parameterization  for  all  of  the  variables  In  Equations  (9) 
and  (10)  mist  be  specified.  The  parameterlzatlons  used  in  the  sporadic-E 
model  which  has  been  developed  can  be  summarized  as  follws: 

a)  Magnetic  field: 

The  magnetic  field  components  were  derived  from  the 
International  Geomagnetic  Reference  Field  (IGRF  1965.0).  Specifically, 
the  three  components  of  the  earth's  magnetic  field  were  computed  in 
geomagnetic  coordinates  from  a series  of  spherical  harmonics  and  their 
derivatives  in  geocentric  coordinates.  The  expressions  used  for 
each  of  the  components  ( x:  east;  y:  north;  z:  vertical)  are  as 
follows: 
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p|J  (cos  9):  Associated  Legendre  function 


a: 


radius  of  the  earth  (6371.2  km) 


h: 


height  above  earth 


9: 


geomagnetic  co latitude 


hn’  V 


IGRF  1965.0  coefficients 


b) 


Ion-neutral  collision  frequency 

The  ion-neutral  collision  frequency,  i.e.,  v^n 
occurring  in  parameter  6^  , was  derived  from  Chapman  (1956).  In  par- 
ticular, Chapman  has  shown  that  v^n  can  be  represented  by  the  follow- 
ing equation: 


vin  = 2.6  x 10'9  M"1/2  N sec'1 


where 


M: 

N: 


ion  mass  (a.m.u.) 

concentration  of  neutrals  (c.g.s.  units) 
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c) 


Neutral  atmosphere: 


An  atmospheric  neutral  density  model  must  be  used  in 
the  process  of  computing  absolute  values  for  the  ion-neutral  collision 
frequency.  Jacchia's  (1971)  neutral  atmospheric  density  model  was 
selected  for  this  purpose.  Jacchia's  particular  parameterization  of 
the  atmosphere  includes  empirical  correlations  between  neutral  den- 
sity and  numerous  geophysical  variables.  The  results  presented  in 
this  chapter  were  computed  using  the  following  values  for  these  geo- 
physical parameters: 


Kp  index: 

Mean  Julian  Day: 

10.7  cm  solar  flux: 

Averaged  solar  flux: 
Solar  declination: 


+ 1 

42495.0  (Vernal  equinox,  1975) 

136  x 10-22  watts  m”2  Hz_1 

155  x 10"22  watts  m“2  Hz“^ 

0° 


c)  Atmospheric  tidal  winds: 

An  important  aspect  of  developing  a numerical  model 
to  study  the  effects  of  tidal  winds  on  sporadic-E  is  the  formulation 
of  a tidal  wind  model.  Classical  solutions  to  the  problem  of  comput- 
ing thermally- forced  tidal  wind  motions  have  been  presented  by 
Wilkes  (1949),  Siebert  (1961),  and  Lindzen  and  Chapman  (1969). 
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However,  these  solutions  do  not  include  the  dissipative  forces  which 
are  particularly  important  in  the  ionospheric  E- region. 


l 


I 


According  to  dissipationless  tidal  theory,  the  amplitude  of 
tidal  winds  increases  exponentially  with  altitude.  Clearly,  this 
is  unrealistic. 

A more  satisfactory  solution  for  the  height  variation  of  the 
tidal  winds  has  been  derived  by  Richmond  (1971)  by  including  the 
dissipative  effects  of  the  Lorentz  force,  molecular  viscosity  and  heat 
conductivity.  The  height  structure  of  the  winds  used  in  the  wind- 
driven  sporadic-E  model  are  those  given  by  Richmond. 

In  the  following  sections  of  this  chapter,  the  theoretical 
temporal  and  spatial  sporadic-E  patterns  generated  by  Richmond's  tidal 
winds  will  be  presented  and  analyzed.  Each  of  the  distinct  tidal 
wind  modes  which  will  be  considered  in  this  study  will  be  designated 
by  two  letters:  n and  m . The  letter  n represents  the  longitudi- 
nal wave  number  (or  the  temporal  periodicity)  of  the  mode  while  the 
letter  m describes  the  latitudinal  structure.  In  particular,  the 
following  standard  conventions  have  been  adopted  for  n and  m in 
this  paper. 


n = 1: 
n = 2: 
m positive: 
m negative: 
m even: 
m odd: 


diurnal  mode 
semidiurnal  mode 
vertically  propagating  mode 
evanescent  mode 

synmetrical  with  respect  to  the  geographic  equator 
asymmetrical  with  respect  to  the  geographic  equator 


1 


! 

' 

■ 
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In  general,  the  various  tidal  wind  modes  which  occur  in  the 
atmosphere  result  from  the  thermal  forcing  produced  by  IR  absorption 
by  water  vapor  in  the  troposphere  and  UV  absorption  by  ozone  in 
the  stratosphere.  However,  the  excitation  and  propagation  of  tides 
is  a complex  problem.  Theoretical  and  observational  studies  indicate 
that  the  (1,1)  and  (1,  -2)  diurnal  modes  and  the  (2,  2)  and 
(2,  4)  semidiurnal  modes  are  the  important  modes  in  the  E-region. 
Therefore,  this  thesis  will  concentrate  on  analyzing  the  spatial  and 
temporal  behavior  of  the  sporadic-E  patterns  produced  by  these 
particular  modes. 

The  vertical  and  latitudinal  structure  for  each  of  the  afore- 
mentioned tidal  wind  modes  is  illustrated  in  Figures  3 and  4 
respectively.  For  the  sake  of  comparison,  the  amplitude  of  each  mode 
has  been  normalized  to  one.  The  phase  profiles  depicted  in  Figure  3 
are  a measure  of  how  the  polarization  of  the  vertically  propagating 
wave  varies  with  altitude.  A graphical  depiction  of  the  east- 
west  and  north-south  components  of  the  important  tidal  wind  modes  is 
presented  in  Appendix  4. 

Using  Richmond's  (1971)  ionospheric  tidal  wind  model,  the 
vertical  ion  velocity  and  vertical  ion  velocity  convergence  patterns 
for  the  (1,  1),  (1,  -2),  (2,  2),  and  (2,  4)  modes  have  been  computed 
for  latitudes  between  10°  N and  70°N.  A summary  of  the  analysis 
of  the  results  for  10°  latitude  increments  along  the  zero  magnetic 
meridian  (=  70°  E longitude)  is  presented  in  Tables  1 through  4. 
Analysis  along  the  zero  magnetic  meridian  was  chosen  to  eliminate 
the  necessity  of  considering  variations  in  the  sporadic-E  patterns 
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Figure  4 Latitudinal  variation  of  Important 
(a-b) 

E- regional  tidal  wind  modes. 


WESTWARD 


GEOGRAPHIC  LATITUDE 


Figure  4(b) 


Semidiurnal  tidal  wind  components 
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Table  1:  Latitudinal  variation  of  the  vertical  Ion  velocity  con- 

vergence produced  by  the  (1,1)  tidal  wind.  (Computa- 
tions based  on  a maximum  N-S  wind  component  of  10  m/s 
at  20°  N.) 


Latitude  (°N) 

10  20 30  40  50  60 


Vertical  propagation 
velocity  of  convergence 

node  (km/hr)  -2.5  -2.5  -2.5  -2.5  -2.5  -2.5 

Maximum  ion  velocity 
convergence. 

(sec" ' x 103)  2.7  4.1  2.9  1.3  .4  .06 

Altitude  of 
maximum  Ion  velocity 

convergence  123  120  123  123  123  119 

(km) 


(See  text  for 
definition.) 


1.72 


1.56 


1.85 


r 


: 


Vertical  propagation 
velocity  of  convergence 
node  (km/hr) 

Maximum  Ion  velocity 
convergence 
(sec_i  x 103) 

Altitude  of 
maximum  ion  velocity 
convergence 
(km) 


(See  text  for 
definition. ) 


Table  3 : 


Latitudinal  variation  of 
vergence  produced  by  the 
tions  based  on  a maximum 
at  200  N.) 


the  vertical  Ion  velocity  con- 
(2,  2)  tidal  wind.  (Computa- 
N-S  wind  component  of  10  m/s 


Latitude  (°N) 


10 

20 

30 

40 

50 

60 

Vertical  propagation 
velocity  of  convergence 
node  (km/hr) 

-10 

-10 

-10 

-10 

-10 

-10 

Maximum  Ion  velocity 
convergence. 

(sec-1  x 103) 

3.2 

5.8 

5.2 

4.2' 

2.5 

1.0 

Altitude  of 
maximum  ion  velocity 
convergence 
(km) 

141 

142 

143 

143 

143 

143 

2.94  3.13 


: 


r 


Table  4:  Latitudinal  variation  of  the  vertical  Ion  velocity  con- 
vergence produced  by  the  (2,  4)  tidal  wind.  (Computa- 
tions based  on  a maximum  N-S  wind  component  of  10  m/s 
at  20°  N.) 


Vertical  propagation 
velocity  of  convergence 
node  (km/hr) 


Maximum  Ion  velocity 
convergence. 

(sec“'  x 103) 


Altitude  of 
maximum  Ion  velocity 
convergence 
(km) 


sn 

(See  text 
for  definition.) 


10 


3.5 


141 


20 


Latitude  (°N) 


1.69 


30 


-5  -5 


40 


1.79 


50 


■5  -5 


60 


3.7  1.1  2.9  2.7  1.3 


142  122  130  137  141 


1.59 


produced  by  magnetic  declination  effects.  To  facilitate  the  compari- 
son of  the  sporadlc-E  patterns,  the  computations  for  each  mode  were 
made  using  a maximum  southward  wind  component  of  10  m/s  at  20°  N. 

In  Tables  1 through  4,  the  parameter  r$n  is  computed  by  dividing 
the  maximum  value  of  shear-produced  ion  convergence  by  the  maximum 
convergence  produced  by  non-shear  mechanisms. 

The  normalized  temporal  and  height  variations  In  the  vertical 
Ion  velocity  produced  by  each  tidal  wind  mode  at  20°  N,  -70°  E are 
presented  in  graphical  form  in  Figures  5(a)  through  5(d).  The  numbers 
in  the  figures  represent  vertical  (upward)  ion  velocities.  The  letters 
represent  vertical  (downward)  ion  velocities.  The  number  5 and  the 
letter  E signify  the  absolute  maximum  normalized  ion  velocities. 

Figures  6(a)  through  6(d)  illustrate  the  normalized  vertical 
ion  velocity  convergence  patterns  corresponding  to  the  data  in 
Figures  5(a)  through  5(d).  The  letter  E in  the  diagrams  represents 
the  maximum  normalized  vertical  ion  velocity  convergence  while  the 
number  5 indicates  the  maximum  ion  velocity  divergence.  Sporadlc-E 
should  form  in  those  areas  where  the  ion  convergence  is  the  greatest. 

The  particular  phases  (i.e.,  local  times)  selected  for  the 
generation  of  Figures  5 and  6 were  specifically  chosen  to  facilitate 
the  comparison  of  these  figures  with  sporadlc-E  observations  presented 
in  the  next  chapter. 

Based  on  the  information  presented  in  Figures  6(a)  through 
6(d)  and  Tables  1 through  4,  the  following  observations  regarding  the 
latitudinal  behavior  of  tidal ly-driven  sporatic-E  can  be  made: 
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Figure  5 
(a-d) 


Normalized  vertical  ion  velocities  gen- 
erated at  20°N,  -70°E  by  important  E- 
reglon  tidal  wind  inodes. 


Legend:  Numbers  represent  vertically  upward  Ion 

motion. 

Letters  represent  vertically  downward  Ion 
motion. 

Note:  Scaling  Is  linear. 
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(U»l)  apinmv 


Figure  5(d)  (2,  4)  mode. 


Figure  6 Normalized  vertical  ion  velocity  con- 
(«-«*> 

vergence/divergence  generated  at  20  N, 

-70°E  by  important  E-region  tidal  wind 

modes. 

Legend:  Numbers  represent  vertical  ion  velocity 

divergence. 

Letters  represent  vertical  ion  velocity 
convergence. 

Dashed  line  designates  an  axis  of  verti- 
cal ion  velocity  divergence. 

Solid  line  designates  an  axis  of  verti- 
cal ion  velocity  convergence. 

Note:  Scaling  is  linear. 
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Figure  6(a) 
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Figure  6(c)  (2,  2)  mode. 
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a)  Sporadic-E  layers  produced  by  the  propagating  tidal 
wind  modes  tend  to  form  in  the  upper  E-region  and  propagate  downwards 
into  the  mid  and  lower  E-region  over  a period  of  6 to  12  hours. 
Therefore,  it  can  be  stated  that  tl daily-driven  E$  layers  generally 
mimic  the  behavior  of  sequentially  sporadic-E. 

b)  The  vertical  propagation  velocity  of  E$  generated 
by  the  evanescent  (1,  -2)  mode  is  a strong  function  of  latitude  (see 
Table  2).  In  particular,  for  latitudes  north  of  30°,  the  numerical 
simulations  Indicate  that  the  (1,  -2)  mode  can  actually  generate 
upward-propagating  sporadic-E  layers. 

c)  The  diurnal  modes  produce  their  greatest  convergence 

at  altitudes  around  125  km.  This  characteristic  is  largely  Independent 
of  latitude. 

d)  In  general,  the  semidiurnal  modes  produce  their  greatest 
convergence  at  altitudes  near  140  km. 

e)  For  the  tidal  wind  modes  under  consideration,  the 
wind-shear  ion  convergence  mechanisms  numerically  dominate  the  non- 
shear convergence  mechanisms  by  factors  ranging  from  1.5  to  3. 

f)  The  overall  convergence  pattern  for  the  (2,  4)  mode 
shifts  6 hours  (90°)  as  30°  latitude  is  traversed.  This  effect  is 
produced  by  the  (2,  4)  tidal  wind  node  situated  near  30°-35°  latitude. 
(See  Figure  4.) 

g)  The  vertical  phase  propagation  velocity  is  distinctly 
different  for  each  of  the  tidal ly-driven  E$  modes.  In  general. 
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sporadic-E  layers  formed  by  the  semidiurnal  tidal  wind  modes  propagate 
vertically  downwards  at  a rate  2 to  4 times  faster  than  E$  layers 
driven  by  the  propagating  (1,1)  diurnal  mode. 

h)  The  magnitudes  of  the  vertical  ion  velocity  conver- 
gences produced  by  representative  atmospheric  tidal  wind  velocities 
( | IJ | ~ 10  to  50  m/s)  are  consistent  with  ion  convergences  actually 
measured  in  sporadic-E  layers  (10  to  10  sec  : MacLeod,  1964). 

i)  The  maximum  wind-driven  ion  convergence  occurs  near 
20°  latitude  for  all  propagating  tidal  wind  modes  which  were  studied. 
This  is  not  true  for  the  evanescent  (1,  -2)  mode. 

Zonal  wind  effects: 

The  effects  of  a vertically-shearless,  east-west  horizontal 

(zonal)  wind  on  the  production  of  sporadic-E  also  has  been  numerically 

evaluated.  The  results  of  the  calculations  for  a vertical  cross- 

section  along  the  zero  magnetic  meridian  (~  -70°  E.  longitude)  are 

shown  In  Figure  7.  This  figure  shows  that  a mean,  west-to-east  zonal 

wind  is  most  effective  in  enhancing  the  formation  of  sporadic-E  in 

the  upper  E-reglon  (~  140  km)  and  at  lower  latitudes.  In  the  same 

figure.  It  can  be  seen  that  zonal  winds  aid  in  the  dissipation  of  E$ 

layers  below  125  km.  The  results  simply  reverse  If  the  effects  of 

a mean  east-to-west  zonal  wind  are  considered.  Peak  convergences 

and  divergences  for  a mean  zonal  wind  speed  of  10  m/s  are  on  the  order 
-3  -1 

of  .2  x 10  sec  at  low  latitudes.  This  Implies  that  zonal  wind 
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Figure  7 Normalized  vertical  ion  velocity  con- 
vergence and  divergence  produced  by 
a mean,  west-to-east,  zonal  wind. 

Legend:  Letters  represent  vertical  ion 

velocity  convergence. 

Numbers  represent  vertical  ion 
velocity  divergence. 


Note: 


Scaling  is  linear. 


(unj)  aprmn 


150. 


» ♦ 

aetieea'iaanB 


l*9i t fHped.ieeBiieB  _ 

»*3.  ccc  iJHBeeeHBbBir 


»*r 

■ 46 

144 

>43 

142 

1 4.1 . 

140 

139 

130 

_11Z. 


136 
135 

134 

Hi 

131 
130 

-m 

127 

120 
125 

124 
123 
122 

121 

<TZ1 

119 

lie 

117 

~IT6 

115 

114 

1 1 1 
110 
IC9 
"IT)  3 

137 
106 

_JLM 

; 04- 

103 

102 

101 

100 


CCCC 
CCcCCC 


eeePBeenBBB 

L‘dOO  JBMJBQH 
BBBBUIWBSBB 
56  3'3H  JffBefl- 


♦ ♦ ♦ 
AAA  AAA AA 4 AAA A 44 A AAA AAA4AA 

A A A A A AAA  AAA  AAA  A A AAA AAA A A_ 

A A A A AAA  AAA  AAA  AAAAAAAAAA 
A A A A A A A A A A A A AAA  A A A A A A 
A AAA AAA AAA AAAAAAAAAA 
AAA AAAAAAAAAA A A A AAA 


fTA  A A AAA  AAAAAAAAAA' 
CCCCCCC:  OB JBBB3BBBB  A A A A A AAA A AA A A AAA A 

CCCCCCCCC  >)BBBOUUBBn  AAA  AAAAAAA  AAAAAA 

_CCCCCCCC 0-3B*»n  IBBHBIJ  A A A AAA  A A A AAAAAA 

CCCCCCCCC  ~ ’ BHUUUH.JOOO  ' A AAA  AAAAAAAAAA 

CCCCCCCCC  UOOOBIIHBBB  A AAAAAAA  A AAA  A 

CCCCCCCCCC  BUHUW1BBBU  AAAAAAAAAAAA 

— t&gftgagt AAAAAAAAAA 


CCCCCCCCCCCCCCCCCCCC 
CCCCCCCCCCCC CCCCCCC 

_CCCCC_C?CCC 

Jo 

UBOOBBB IB 

AAA  BBUeH'inoaU  I'lnilBHB'IBBHOtJOOUBBO 

A A AAA  BnH'IHUOHBOUHUO 

— rTTT 


OUBOBOOtlUO 
BHOBO  1HIM1B 
'oionnoHiioo 
BuliOeOBBBBQ 
UBBOllOHOHiilittB 


A A A A A A A AVaA 
A A AAA AAAAAA 
AAAAAAAAAA 
AAAAAAAAAA 
AAAAAAAAAA- 
AAAAAAAAAA 
AAAAAAAAAAA 
AAAAAAAAAAAA 
mAAAAAAAAAAX* 


flnkXSXA 

222  111  11  AAAAAAAAAAAA  A AAA A A A A AAA A AAA AAA AA 

3 2222  Hill  AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

331  2232  1 1 111  11  AA  AA  A AAA A A A A AAA A AAA 

T44-  “ 2T3T3  JT225Z 1 nm  fl 

4444  3333  272222  1111111111 

556  4444  3333  222222  1111111111 

5555  44444  13133  227222  1111111111 

— 5=353 


665656 

5555555 
555555 
"5b 


-rrrr 

44444 
44444 
4444444 
“AAA  A 4*6 — 


44464444 
444444446464 
4444 


TT7TT 
333J3J 
33  333 
333333 
TI333T 


3333333 
3333333 
33333133 


m7T?. n r i nn  1 1 

222222  111111111 

22222  lllltntt 

222222  lllllltl 

223222 mrrrrr 

222222  llllllll 


222222 

2222222 


i u mi  i 
l ilium 


rrrnrrrrTri 22222222 riTrriTn" 

333J333J33333J33  222222222  1111111111 

3 33333  22272222222  1111111111 

2222222722222  lllllllltlll 

“222 ’277’7  272777727777277 II  I ITlTHTni 

22222222222222232  11111111111111 

111  III  till  II  till  I 

1111111  1 1 11  1 1 I till  1 1 

1 1 1 1 m u i m 1 1 nT  mum 
ii  i un  1 1 1 ii  i it  m 1 1 mi  i m mu  i it 
i n in  in  i u i in  1 1 mu  1 1 in  i ui 
i u i i in  u u in  i u i mi  m 

~ it  i vn  mi  i n mi 

6 » 4 4 4 ♦ 

10  2 9 30  40  SO  60 


70 


Latitude  (°N) 


effects  are  approximately  an  order  of  magntltude  smaller  than  the  Ion 
convergences  generated  by  the  propagating  tidal  winds. 


S unwary 

In  this  chapter,  the  details  of  the  construction  of  a 
tl daily-driven,  numerical  sporadlc-E  model  were  discussed,  and  results 
of  several  simulations  were  presented.  The  results  indicate  that 
E$  layers  theoretically  can  be  produced  by  vertical  ion  convergences 
generated  by  atmospheric  tidal  winds.  Furthermore,  the  simulations 
Indicate  that  each  different  tidal  wind  mode  produces  a sporadic-E 
layer  with  a distinct  set  of  behavioral  characteristics.  Therefore, 
by  studying  observed  spatial  and  temporal  patterns  of  sporadic-E,  It 
may  be  possible  to  Identify  and  study  the  Important  diurnal  and  semi- 
diurnal atmospheric  tidal  wind  modes.  This  particular  problem  will 
be  considered  In  the  next  chapter. 
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IV.  COMPARISON  OF  THEORETICAL  RESULTS  WITH  SPORADIC-E  DATA 


In  the  previous  chapter,  the  spatial  and  temporal  ion  con- 
vergence and  divergence  patterns  for  selected  atmospheric  tidal  wind 
modes  were  presented.  In  this  chapter,  the  results  of  these  theoreti- 
cal numerical  computations  will  be  compared  to  published  sporadic-E 
data. 

The  task  of  evaluating  theoretical  E$  model  predictions  is 
greatly  hampered  by  the  paucity  of  useful  sporadic-E  data.  In  general, 
most  published  sporadic-E  data  are  climatological  in  nature.  Hence, 
detailed  evidence  which  might  confirm  specific  features  predicted  by  a 
sporadic-E  model  is  usually  not  available.  The  evaluation  of  E$ 
model  results  is  further  limited  by  current  ionosonde  data  reporting 
standards.  Under  the  present  URSI  ionogram  reporting  system,  the 
virtual  height  of  only  one  sporadic-E  layer  can  be  reported  during 
any  given  observation  period  (usually  1 hour).  However,  as  noted  in 
the  first  chapter,  the  presence  of  more  than  one  E$  layer  in  the 
ionosphere  is  coimton  (see  Figure  1,  profile  1).  Since  sporadic-E 
layers  in  the  lower  E-region  tend  to  be  more  concentrated  in  compari- 
son to  Es  layers  in  the  early  stages  of  formation  in  the  upper  E- 
region,  the  lower  layers  are  reported  and  the  higher  layers  are 
ignored.  Also,  since  sporadic-E  layers  with  critical  frequencies 
less  than  .5  MHz  are  not  reported  under  the  URSI  system,  the 
embryonic  behavior  of  high-level  sequential  sporadic-E  is  largely 
undocumented.  The  usefulness  of  Central  Radio  Propagation  Laboratory 
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(CRPL)  sporadic-E  data  Is  also  severely  limited  by  the  poor  virtual 
height  resolution  (5  km)  used  In  the  reporting  system.  This  coarse 
height  resolution  particularly  Inhibits  the  accurate  computation  of 
the  vertical  phase  propagation  velocity  of  sequential  sporadlc-E. 
Finally,  since  most  lonosonde  stations  have  been  established  for  the 
primary  purpose  of  observing  and  reporting  the  behavior  of  photo- 
chemically-produced  E-  and  F-reglon  Ionization,  It  Is  not  uncotimon 
to  find  many  stations  operating  on  daytime  schedules  only.  Therefore, 
comprehensive  studies  of  night-time  sporadlc-E  are  not  easily  conducted. 

In  order  to  circumvent  the  limitations  of  the  present  world- 
wide lonosonde  data  reporting  system,  MacDougall  (1974a, b)  has  hand- 
analyzed  numerous  unreduced  Ionograms  for  the  expressed  purpose  of 
accurately  determining  the  virtual  heights  of  sporadic-E  layers  as  a 
function  of  time.  The  dots  in  Figure  8 represent  occurrences  of 
sporadic-E  over  Jamaica  (20°  N,  -70°E)  as  determined  by  MacDougall. 

The  lines  in  Figure  8 have  been  transcribed  from  Figures  6(a)  through 
6(d).  They  represent  the  theoretically- computed  axes  (or  planes)  of 
vertical  ion  velocity  convergence  which  are  formed  by  using  Richmond's 
(1971)  tidal  winds  as  the  forcing  function  in  the  numerical  sporadic-E 
model  presented  In  the  previous  chapter. 

At  this  point,  It  is  important  to  recognize  that  the  analysis 
presented  In  Figure  8 represents  a comparison  between  the 


Figure  8 Analysis  of  sporadic-E  virtual  height 
data  for  Jamaica  (20°N,  -70°E)  for  1-15 
March  1965. 
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virtual  heights  of  sporadic-E  layers  and  theoretically-derived, 
geometric  heights  of  maximum  ion  convergence.  In  general,  they  are 
not  equal.  However,  considering  the  relatively  low  total  electron 
content  of  the  lower  E-region,  it  is  reasonable  to  approximate  virtual 
heights  by  actual  heights  in  the  following  analysis. 

The  data  presented  in  Figure  8 suggest  that  there  is  a strong 
connection  between  atmospheric  tidal  winds  and  the  occurrence  and 
behavior  of  sporadic-E.  In  particular,  the  following  features  should 
be  specifically  noted: 

a)  Certain  virtual  height  patterns  of  sporadic-E  can  be 
specifically  attributed  to  ion  convergence  generated  by  individual 
tidal  wind  modes. 

b)  The  phases  of  the  individual  tidal  wind  modes  remain 
relatively  constant  during  the  two-week  equinox  period  depicted  in 
Figure  8. 

c)  Sporadic-E  at  low  latitudes  is  largely  sequential  in 
nature,  i.e.,  the  layers  descend  from  the  upper  to  the  middle  and  lower 
E-region  over  a period  of  6 to  12  hours.  Furthermore,  the  rate  of 
descent  is  distinctly  different  for  each  tidal  wind  mode. 

d)  The  frequency  of  occurrence  of  E$  is  highest  where 
the  ion  convergence  generated  by  two  or  more  modes  is  mutually  rein- 
forcing, e.g.,  16  to  18  LT  at  130  km  in  Figure  8. 

e)  Metallic  ions  are  apparently  lost  from  the  E-region  as 
they  diffuse  into  the  mesosphere  through  the  turbopause  (*  100  km). 
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The  lines  in  Figure  8 represent  only  the  axes  of  maximum  ver- 
tical Ion  velocity  convergence  for  the  selected  tidal  wind  modes. 
However,  corresponding  to  each  axis  of  Ion  convergence  Is  an  axis  of 
vertical  Ion  velocity  divergence  (see  Figures  6(a)  through  6(d)). 

The  lines  of  maximum  Ion  divergence  corresponding  to  the  lines  of 
convergence  drawn  In  Figure  8 are  shown  separately  In  Figure  9. 

The  axes  of  Ion  divergence  shown  In  Figure  9 explain  several  Interest- 
ing features  of  Figure  8. 

a)  The  low  frequency  of  occurrence  of  sporadic-E  during 
early  morning  hours  (0-10  LT)  can  be  partly  explained  by  the  mutual 
reinforcement  of  the  combined  divergence  patterns  of  t e four  tidal 
wind  modes.  The  fact  that  much  of  the  neutral  metal  ^bris 
entering  the  atmosphere  during  the  night  is  not  photc  , ed  until 
sunrise  is  undoubtably  an  additional  contributing  factor. 

b)  The  "ion  notch"  (i.e.,  depression  of  E$  virtual 
height)  seen  at  120  km  between  7 and  12  LT  in  Figures  8 and  10  can 

be  conveniently  explained  by  the  axis  of  Ion  divergence  for  the  (1,1) 
mode  (see  Figure  9). 

Figure  10  shows  the  occurrences  of  sporadlc-E  observed  during 
the  second  half  of  March,  1965  at  Jamaica  By  carefully  comparing  the 
data  In  Figures  8 and  10,  several  noteworthy  observations  regarding 
the  temporal  behavior  of  sporadic-E  can  be  made: 

a)  The  phases  of  all  of  the  vertically-propagating  modes 
((1,  1),  (2,  2),  (2,  4))  has  shifted  approximately  1 hour  during  the 
one  month  period  of  observation.  This  suggests  that  the  phases  of 
the  tidal  winds  mey  be  functions  of  the  time  of  year. 
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Figure  9 Axes  (or  planes)  of  Ion  divergence  cor- 
responding to  the  axes  of  ion  convergence 
depicted  in  Figure  8. 
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Figure  10  Analysis  of  sporadic-E  virtual  height  data 
for  Jamaica  (20°N,  -70°E)  for  15-31  March 
1965. 


b)  Production  of  sporadic-E  by  the  (2,  4)  tidal  wind  mode 
declined  during  the  second  half  of  March  1965.  This  implies  that  the 
amplitudes  of  the  tidal  winds  are  not  constant  with  respect  to  time. 

c)  There  is  little  apparent  change  in  the  frequency  of 
occurrence  of  sporadic-E  generated  by  the  diurnal  (1,1)  and  (1,  -2) 
wi nds . 

Assuming  that  the  phases  for  the  prominent  equinox  tidal  wind 
modes  have  been  correctly  determined  by  the  analysis  of  MacDougall's 
Jamacian  data  (Figures  8 and  10),  it  should  be  possible  to  further 
test  the  validity  of  the  sporadic-E  numerical  model  computations 
by  attempting  to  extend  the  model's  predictions  to  middle  and  high 
latitude  stations.  This  has  been  done.  In  particular,  sporadic-E 


virtual  height  patterns  for  various  Far  East  ionosonde  stations  have 


been  analyzed.  The  equinox  sporadic-E  data  for  Akita,  Japan 
(39°N,  140°E)  is  presented  in  Figure  11.  Although  the  data  lacks  the 
resolution  and  detail  present  in  MacDougall's  data,  it  still  illustrates 
several  important  features  of  mid- latitude  sporadic-E: 

a)  The  sporadic-E  virtual  height  patterns  seen  at  mid- 
latitudes (30°  - 60°)  are  not  as  distinct  as  those  produced  at  low 
latitudes  (10°  - 20°).  This  is  partly  due  to  the  fact  that  tidally- 


drlven  ion  convergence  reaches  a maximum  near  20°  latitude  for  the 
important,  non-evanescent.  E-region  tidal  wind  modes  (reference  Tables 
1 through  4). 

b)  Due  to  the  presence  of  the  (2,  4)  wind  nodes  near 
* 30°  latitude  (see  Figure  4),  the  (2,  2)  and  (2,  4)  Ion  convergence 
and  divergence  patterns  are  out  of  phase  with  respect  to  each  other 
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Figure  11  Analysis  of  sporadlc-E  virtual  height 
data  for  Akita,  Japan  (39°N,  140°E)  for 
March  1974. 


Figure  12  Axes  (or  planes)  of  Ion  divergence 
corresponding  to  the  axes  of  Ion 
convergence  depicted  In  Figure  11. 
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at  latitudes  greater  than  30°.  Consequently,  they  destructively 
Interfere  with  one  another  (see  Figures  11  and  12).  Therefore,  the 
mid- latitude  ion  convergence  patterns  produced  by  the  semidiurnal 
tides  are  not  as  distinct  as  the  patterns  depicted  in  low-latitude 
E$  data,  e.g. , Figures  8 and  10.  (*See  footnote  at  end  of  chapter.) 

c)  The  (1,-2)  mode  produces  ascending  E$  layers  at 
mid-latitudes.  This  behavior  is  contrary  to  the  descending  ion  con- 
vergence patterns  produced  by  the  other  tidal  wind  modes.  As  such, 
the  (1,  -2)  ion  convergence  pattern  adds  considerable  complexity  to 
depictions  of  sporadic- E occurrences  at  middle  and  high  latitudes. 

d)  The  ion  convergence  patterns  produced  by  the  (1,1) 
and  (2,  4)  modes  at  mid-latitudes  are  nearly  in  phase  during  the 
morning  hours  (6-12  LT).  Accordingly,  the  strong  convergence  gener- 
ated by  the  mutual  interaction  of  these  modes  is  clearly  evident  in 
the  sporadic-E  virtual  height  patterns  for  Akita  (Figure  11). 

As  an  aside,  it  should  again  be  reemphasized  that  the  phases 
of  the  tidal  winds  used  in  the  comparison  of  the  theoretical  and 
observational  results  presented  in  this  chapter  were  empirically 
chosen  to  maximize  the  agreement  between  MacDougall's  data  and  the 
E$  numerical  model  predictions.  Nonetheless,  the  semidiurnal  tidal 
wind  phases  used  for  the  analysis  of  the  sporadic-E  patterns  in 
Figures  8 and  10  are  in  good  agreement  with  wind  phases  independently 
deduced  from  the  analysis  of  magnetic  field  variations  produced  by 
ti daily-driven  ionospheric  currents.  (Unpublished  work  at  UCLA.) 
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Specifically  the  (2,  4)  phases  are  virtually  equal  while  the 
(2,  2)  phases  differ  by  1.5  hours  (*  22°).  The  determination  of 
the  phase  of  the  (1,  -2)  tide  from  sporadic-E  data  is  somewhat 
more  uncertain.  Although  the  axis  of  maximum  ion  convergence  for  the 
(1,  -2)  mode  appears  to  be  properly  placed  in  the  Jamacian  data,  it 
seems  that  a better  agreement  with  the  Akita,  Japan  data  can  be  made 
If  the  phase  is  shifted  minus  4 hours.  The  overall  difficulty  In 
accurately  determining  the  phase  of  the  (1,  -2)  tidal  wind  is  in  part 
due  to  the  fact  that  the  (1,  -2)  mode  is  not  particularly  effective 
in  producing  ion  convergence  in  the  E-region  (see  Table  2).  Hope- 
fully, this  problem  can  be  overcome  through  the  analysis  of  more  mid- 
latitude ionosonde  data  since  the  (1,  -2)  mode  produces  its  maximum 
ion  convergence  near  40°  latitude. 

Summary : 

The  purpose  of  this  chapter  was  to  clearly  demonstrate  the 
fact  that  atmospheric  tidal  winds  play  a distinct  and  important  role 
in  the  formation  and  vertical  propagation  of  sequential  sporadic-E. 

In  particular,  it  was  shown  that  Richmond's  (1971)  tidal  wind  formula- 
tion, when  used  as  the  forcing  function  in  a numerical  sporadic-E 
model,  correctly  predicts  the  broad  features  of  equinox  E$  for 
two  widely-separated  stations  situated  at  different  latitudes. 


The  statistical  significance  of  the  Akita,  Japan  Es  data  is 
presently  being  improved  by  including  Vernal  equinox  data  for  several 
years. 
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V. 


CONCLUSIONS 


The  most  interesting  aspects  of  the  mathematical  analysis  and 
observational  evidence  presented  in  the  previous  chapters  can  be  briefly 
summarized  as  follows: 

a)  The  gross  features  of  the  convergence  of  metallic  ions 
in  the  Ionospheric  E-region  (sporadic-E)  can  be  adequately  explained 
by  considering  a steady-state  equilibrium  between  neutral  wind  motions 
and  the  resultant  collisional  and  geomagnetic  forces  experienced  by 
the  ions  In  the  neutral  gas. 

b)  Representative  tidal  wind  velocities  generate  sufficient 
ion  convergence  to  account  for  the  observed  formation  of  sporadic-E 
layers. 

c)  Each  tidal  wind  mode  generates  a theoretically-distinct 
pattern  of  vertical  ion  velocity  convergence.  Each  characteristic 
pattern  is  a function  of  geographic  and  geomagnetic  latitude. 

d)  Occurrences  of  sporadic-E  attributable  to  ion  conver- 
gence generated  by  individual  tidal  wind  modes  have  been  Identified  in 
composite,  low-latitude  E$  data.  As  such,  it  has  been  possible  to 
estimate  the  phases  of  the  important  E-region  tidal  wind  modes. 

e)  Having  estimated  the  phases  and  approximate  amplitudes 
of  the  tidal  winds  at  a given  location.  It  is  theoretically  possible 
to  describe  the  gross  temporal  patterns  of  sporadic-E  which  should  be 
observed  at  other  latitudes  and  longitudes.  This  hypothesis  has  been 
put  to  a test  through  the  analysis  of  equinox  sporadic-E  data  from 
Jamaica  and  Japan. 
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At  present,  the  predictive  capabilities  of  the  numerical 
sporadic-E  model  which  has  been  presented  seem  promising.  However, 
the  ultimate  usefulness  of  the  model  will  have  to  be  determined  through 
the  careful  analysis  of  more,  global  sporadic-E  data. 
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Derivation  of  the  ion  velocity  resulting  from  mech 
ical  and  electrostatic  forcing  of  a collisional, 
magnetic,  weakly- ionized  plasma  (e.g.,  the 
ionospheric  plasma). 


In  general,  the  equation  of  motion  for  ions  in  the  earth's 
ionospheric  E-region  can  be  written  as 


raini  (ft-  = miVin 


(U  - Vi)  + 


±i  x i \ 

~~Z )+  minifl-  “ VPi 


(1.1) 


where:  m^ 

= 

ion  mass 

ni 

= 

ion  concentration 

h 

= 

ion  velocity 

vin 

= 

ion  neutral  collision  frequency 

U 

= 

neutral  wind  velocity 

= 

ion  charge 

E 

= 

electrostatic  field  vector 

B 

= 

magnetic  field  vector 

C 

= 

speed  of  light 

1 

= 

gravitational  vector 

VP1 

= 

ion  partial  pressure  gradient 

For  externally- forced  ion  motions  In  the  ionospheric  E-region, 
the  Lorentz  and  collisional  forces  dominate  the  gratitational  and 


ion  partial  pressure  gradient  forces.  (See  main  text  for  a detailed 
discussion.)  Therefore,  by  neglecting  the  Ion  partial  pressure 
gradient  and  gravitational  terms,  the  steady-state  equation  of  motion 
simplifies  to: 


/ L x B \ 

Vin  W-  ill  * Vi  (l*  -V- ) 


n^n.-v.. 


which  can  be  rewritten  as: 


q,|B|  /CE  V.  x B\ 

vin (^-  V + Tqr  ynj  + nj-7 


(1.2) 


By  defining  the  following  quantities: 


e|B| 

<*>i  = : i°n  gyro  frequency 


e = 


CE 

IF 


: reduced  electric  field 


b = 


B 

nr 


unit  vector  parallel  to 
magnetic  field  vector 


Equation  (1.2)  becomes: 


in 


— (U  - ) + e + (V^  x b)  « 


(1.3) 
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Define:  B,  = ^ . 

1 U),. 

Using  this.  Equation  (1.3)  can  be  rewritten: 

Bl  (H  - Jtj)  + e + (V,,  x b)  - 0 
Crossing  Equation  (1.4)  with  b , we  have: 

3i  (U  x b)  - 0.  (V.  x b)  + e x b + (V^.  x b)  x b = 0 (1.5) 

Next,  recall  that: 

i 

(Vj  x b)  x b = - Vf  (b-b)  + (V. -b)  b 

A /V 

or,  since  b*b  = 1 : 

(Ij  x b)  x b = - Vf  + (V.-b)  b (1.6) 

Substituting  Equation  (1.6)  into  Equation  (1.5): 

6i  (U  x b)  - 0i  (V.  x b)  + e x b = ^ - (V^b)  b (1.7) 

From  Equation  (1.4),  we  note  that: 


77 


±i  - j:  (e-b)b  - (VixsSN^b-  (U-b)  b 


0.10) 


Equation  (1.10)  now  can  can  be  simply  solved  for  V . , the  vertical  ion 
velocity  for  a wind-driven,  collisional,  magnetic  plasma: 


li 


+ ei  (U  x b)  + (U-b) 


A A I 

U-b)  bl 


[*?][ 


Bj£  + e x b + T-  (£*b)  b 


J 


(1.11) 
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Note  that  the  r.h.s.  of  Equation  (1.11)  has  six  terms.  The  first 
three  represent  the  vertical  ion  velocity  contribution  from  neutral 

wind  forcing.  The  second  three  terms  represent  the  contribution  from 
electric  fields. 


APPENDIX  2 


Derivation  of  the  Ion  velocity  divergence/convergence 
resulting  from  mechanical  forcing  of  the  earth's 
ionospheric  plasma. 


The  ion  velocity,  V.  , attributable  to  the  collisional 
coupling  between  neutral  wind  motions  and  ions  in  a magnetic,  weakly- 
ionized  plasma  can  be  readily  Identified  as  the  first  three  terms  in 
Equation  (1.11),  Appendix  1: 


u pi  pi  - i 

-i  = , . „2  ^ + T 7 (U  x b)  + — 1 


1 + *J  L 1 + 0 


(U*b)  b (2.1) 


Taking  the  divergence  of  Equation  (2.1): 


0- 

V-V.  = V — , 
1 + 6' 


0-  6- 
* V.  + 7 (V-U)  + v — ■ — • (u  X b) 

[i  ♦ »?J  77^- 


* 771?  '•(!!  x b)  ♦ 7 — (u-b)  -b  + — !—y  (u-b)-  Vr 

L1  ei  [i  + ef  [ l +ej  " J ^ 


Expanding  terms  using  vector  Identities: 


Substituting  Equation  (2.3)  Into  Equation  (2.2),  « havo  the  expression 
for  the  ion  velocity  divergence/convergence  arising  from  ion  motions 
collisionally  Induced  by  neutral  wind  motions: 
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b • c (b • v)uj 


+ (u-b)  V 


T + ft? 
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APPENDIX  3 


Derivation  of  the  ion  velocity  divergence/convergence 
resulting  from  the  electrostatic  forcing  of  the 
earth's  ionospheric  plasma. 


The  ion  velocity,  , attributable  to  motion  induced  by  an 
electrostatic  field  in  a collisional,  magnetic,  weakly-ionized  plasma 
can  be  readily  identified  as  the  last  three  terms  in  Equation  (1.1), 
Appendix  1 . 


\ 


Expanding  terms  using  vector  identities: 


V 


e + 


Pi 

1 + 0? 


(V-e)  + 7 


Expanding  the  last  term  in  Equation  (3.2): 


(V(e-b))-b 


Substituting  (3.3)  back  Into  Equation  (3.2),  we  have  the  expression 
for  the  Ion  velocity  divergence/convergence  arising  from  Ion  motions 
Induced  by  electrostatic  fields: 


APPENDIX  4 


Graphical  depiction  of  the  vertical  and  temporal 
structure  of  the  N-S  and  E-W  components  of 
the  important  E-region  tidal  wind  modes. 

Legends:  (a)  For  E-W  winds: 

Numbers  represent  eastward  winds. 
Letters  represent  westward  winds. 

(b)  For  N-S  winds: 

Numbers  represent  northward  winds. 
Letters  represent  southward  winds. 


Note : 


The  following  diagrams  were  generated 
with  the  same  phases  which  were  used 
in  computing  the  Ion  velocity  and  con- 
vergence/divergence diagrams  presented 
in  the  text.  The  effects  corresponding 
to  phases  other  than  those  shown  here 
can  be  computed  by  shifting  the  diagrams 
minus  1 hour  for  each  +30°  phase  shift 
in  the  semidiurnal  winds  and  minus  1 
hour  for  each  +15°  phase  shift  In  the 
diurnal  winds. 
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APPENDIX  5 


Scale  analysis  of  vertical  ion  velocity  convergence/divergence 
produced  by  representative  tidal  winds  at  40°N,  -70°E. 

A scale  analysis  of  wind-  and  electrostatically-driven  verti- 
cal ion  velocity  convergence  is  greatly  simplified  if  the  analysis  is 

A 

performed  at  a location  along  the  zero  magnetic  meridian  where  bx 
equals  zero.  By  including  the  simplifying  assumptions  presented  in 
Chapter  II,  the  vertical  ion  velocity  convergence  equation  for  the 
zero  magnetic  meridian  becomes: 


a-z  _ _b  l : b 

32  62  y ' 8z  * 2 


where 


✓N.  ^ Q % 

+ -x  by  fz  ■ ^y  by  bz  az 


9 1 


r h 9 1 


B = v.n/a). 


U : neutral  wind  field 


V : ion  velocity 


b : unit  magnetic  vector 


For  the  purpose  of  scale  analysis,  a constant,  horizontal, 
height-independent,  E-W  electric  field  of  1 mv/m  has  been  specified. 
This  equates  to 


The  scale  analysis  has  been  performed  for  the  five  important 
ion  convergence  mechanisms  for  the  altitudes  of  120  and  150  km  using 

l^max^  **  50  m/s.  The  results  follow. 


(a)  Scale  analysis  for  120  km  (altitude)  tidal  winds: 
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